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Abstract. We investigate the density of states for a photon gas confined in a 
nonmagnetic metamaterial medium in which some components of the permittivity 
tensor are negative. We study the effect of the resulting hyperbolic dispersion relations 
on the black body spectral density. We show that for both of the possible wave-vector 
space topologies, the spectral density vanishes at a certain frequency. We obtain the 
partition function and derive some thermodynamical quantities of the system. To 
leading order, the results resemble those of a one or two dimensional photon gas with 
enhanced density of states. 
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1. Introduction 

The density of states is a basic factor governing the physical properties of many 
statistical systems. This can be seen, for instance, by considering the effect of electronic 
density of states on conductivity or other properties of condensed matter systems. 
Similarly, there are many properties of optical systems which depend on photonic density 
of states. This fact has been used in designing photonic devices with desired properties 
by altering their photonic density of states through a variety of techniques [TJ [2j [3j H] , 
and more recently, by using metamaterials [5] . The later method relies on the properties 
of metamaterials, which, roughly speaking, are media with negative indices of refraction 
in certain frequency ranges. 

The propagation of electromagnetic waves in media with negative permittivity and 
or permeability was studied decades ago by Veselago [6] . This results, depending on the 
signs of the permittivity and the permeability, in the emergence of so-called evanescent 
or backward-propagating electromagnetic waves (for a brief review, see e.g., Ref. [7] 
and for a more extensive review see Ref. [8]). Following the experimental realization of 
such media [9j QUI El 1121 13], the subject has attracted a lot of interest in both science 
and technology. One can mention in this regard, on the practical side, the construction 
of various microwave devices ranging from antennas to band-pass filters and lenses, see 
e.g., Ref. [H] and references therein. On the other hand, from the more formal point of 
view, metamaterials are considered as arenas in which one can mimic certain theoretical 
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models such as dynamic spacetime, string theory D-branes and noncommutativity [15J, 
spacetimes with a metric signature transition [16J, Schwarzschild (anti-) de Sitter black 
holes [TF], and motion around celestial bodies p]S] . 

The present work aims to study the thermodynamics of a photon gas confined within 
a metamaterial medium, which usually obeys hyperbolic dispersion relations [19]. This 
is partly motivated by the recent interest in relativistic gases with modified dispersion 
relations. Such modifications have been considered within several contexts, namely, in 
Refs. [201 ED E2J E2J EI], in connection with Lorentz invariance violating models within 
the context of quantum field theory; in Ref. [25], in the context of quantum gravity; in 
Ref. [26] in the context of high energy astrophysics, in Refs. [271 EH] in the framework 
of noncommutative field theories, in Ref. [29], in the context of models with minimal 
length; and in Ref. [30], in the framework of the Hofava gravity. Another motivation 
behind this work is the interest in developing different probes into the properties of 
metamaterials. As an example, spontaneous emission near hyperbolic metamaterials 
has been considered in Refs. [3U 132] . 

The hyperbolic dispersion relation associated with nanostructured metamaterials 
and its application in engineering photonic density of states has been discussed in Ref. 
[5]. Here, we study the statistical mechanics of a photon gas with such dispersion 
relations for different wave-vector space topologies and obtain several new results. In 
particular, we show that they exhibit interesting properties such as the vanishing of 
density of states at certain frequencies and behaving like lower dimensional usual photon 
gases but with enhanced density of states. 

In the next sections, after a brief review of the dispersion relation for the 
propagation of electromagnetic fields in an anisotropic medium, we show how hyperbolic 
dispersion relations could arise as a result of negative permittivities. We then obtain an 
expression for the photonic density of states in a rather simple model. We apply this 
to study the black body radiation and show that this results in novel properties such 
as exhibiting the behaviour of a photon gas in lower dimensional ordinary media. We 
study the thermodynamics of such media by calculating the partition function in the 
grand canonical ensemble. We conclude by discussing the results. 

2. Hyperbolic dispersion relations 

The propagation of electromagnetic waves is governed by the well-known Maxwell 
equations 
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in which D = e-E and B = /i-H. In general, neglecting medium losses, the permittivity e 
and the permeability /i are symmetric tensors. Here, we consider an anisotropic dielectric 
medium with 



(5) 



where e XjV)Z are all constant, and seek plane wave solutions of the above equations in the 
absence of the electric charge p and the current J. Thus, we set 

E(r,t) = E e ik - r ^ (fc) ' (6) 

where the real part is to be taken. Inserting this back into the Maxwell equations (J2J) 
and (j3J and making use of Eq. (pQ), one can show that the following relation holds 

k 2 k,, k 



+ 



k 2 



k 2 - 
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(7) 



in which 



x, y, z), and c 



-^= is the speed of light in vacuum. If we 

confine ourselves to the special case of a uniaxial medium with e a 
the following simpler relation 



e y , we can obtain 



k\ + k 2 



+ — 



For ordinary media, e x and e z are positive and can be constant. However, for 
metamaterials, either of these can have negative values. Also, in general, metamaterials 
are always highly dispersive media, i.e., they exhibit both temporal (frequency 
dependence) and spatial (wave- vector dependence) dispersions [33J. Thus, the explicit 
forms of the permittivity components depend on the structure of the medium under 
consideration. Since the spatial dispersion effects are usually much weaker than the 
temporal ones, here, for convenience, we consider a dissipationless medium and ignore 
the spatial dispersion. In the presence of dissipation, which is the case for the metal- 
dielectric metamaterials, the permittivity components are complex and the above simple 
hyperbolic dispersion get replaced by hyperbolic-like dispersion relations [53] . However, 
we expect the results obtained here to be qualitatively valid for low-loss media as well. 
We consider a medium with 



(8) 
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(9) 
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in which Q is some cutoff frequency. With this choice, we obtain negative values for e x 
for frequencies under the cutoff value. This characterizes a nonmagnetic metamaterial 
with negative xx and yy components of the permittivity. Thus, Eq. ([8]) reduces to 



J-2 , u2 _ U 



k: 
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An example of such dispersion relations has been used in Ref. [32] (see also Ref. [36J) 
to study the physics of vacuum in a strong magnetic field in which the permittivity 
components are 

1 + a 

E-xx ^yy j ~ 

and 

where u s is a constant frequency proportional to the inverse London penetration depth 
and < a < 1 is also a constant. 

In this work, we confine ourselves to the case of nonmagnetic metamaterials. To 
take metamaterials with a magnetic response into account, one would be faced with a 
more complicated picture, in which a variety of mode-dependent dispersion relations 
(classified in Ref. |37]) should be considered. 

3. The density of states 

To study the statistical mechanics of the system, we first need to calculate the density 
of states. This can be obtained, as usual, from the wave- vector space volume of the shell 
enclosed by surfaces defined by u and u + du. Obviously, this volume is infinite, giving 
rise to a divergent density of states. However, in practice, there is a bound on the values 
of k which is imposed by the physical properties of the medium. An example of this is 
discussed in Ref. [5] in which the metamaterial under study is fabricated via a special 
nanopatterning with a characteristic patterning scale a. This puts an upper cutoff value 
on k which is of the order of — . Thus, we need to calculate the above volume taking 
the constraint k < k c into account, k c being the cutoff value. The volume inside the 
hyperbolic surface described by Eq. (ITT1) (depicted in Fig. [T]) and the planes k z = ±k , 
where ko = , equals 

where we have used k c c = Q. From this, we obtain 

Basically, this gives the density of states in the interval [oj,u + du]. However, 
there are some exotic properties associated with this relation, namely, the appearance 
of negative values. This is because of the presence of a fixed cutoff value for k. In certain 
regions, by increasing u, the value of k Q decreases so the net effect is a decrease in the 
enclosed volume. The behaviour of ^ is depicted in Fig. [2j 

The number of states in the interval between ui and oj + du may be obtained by 
multiplying a factor of by the absolute value of ^£duj. Here, V is the volume of 
the medium and the factor | is inserted to account for two polarization states of each 
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Figure 1. The surface spanned by the dispersion relation k 2 x + ky — q2_ u 2 k-1 = in 
the wave- vector space. The density of states is given (up to a proportionality constant) 
by the volume enclosed by the hyperboloid and the planes k z = ±fco- 




0.8 1 



Figure 2. The function v' = 4^ (in units of %^ ) in terms of ui (in units of fi). 



photon, and to include the positive octant only. Thus, the density of states p{u) is given 
by 

V 



p(u)du 



7 co 



AO? 



lOfiV + 3^\duj. 



(12) 



3n 2 c 3 W 

For oj > Q, the medium obeys the usual spherical dispersion relation, and we have 
p(uj)du> = -^-^u 2 duj. 

An immediate effect of the above modified relation for the density of states can 
be seen by studying the black body radiation. To see this, we consider a metamaterial 
medium of the type described above and a gas of photons in equilibrium with the 
medium at a temperature T = — ^ (kg being the Boltzmann constant). The energy 
density of the radiation emitted inside the medium at this temperature is given by the 
following expression 

1 hwp(u)du 



u(uj)dt>j 



y e /3hui _ i 



(13) 
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where u(u) is the spectral density. For the usual spherical dispersion relation, this 
leads to the well-known Stefan-Boltzmann law. However, using the hyperbolic density 
of states, Eq. ( TT2l) . the result would be different. In practice, we need not to use the 
whole expression Eq. (fl2j) here, but rather a simplified version would be sufficient. To 
see this, we first rewrite the above expression in the following form 



u(x) 



\AX A -WX 2 x 2 + 3x 4 | 



x 



(14) 



37r 2 c 3 ft 3 /3 4 X 3 e* - 1 

where x = f3hu and X = (3h~Q. Now, we note that for a medium with nanoscale 
characteristic patterning length a ~ 10~ 8 m, for room temperatures we have X ~ 10 3 , 
which means that we can practically confine ourselves to the region i<I. Thus, we 
obtain 



AX 



U{X) = 3 7rW/J*e*-l ( 15 ) 
The functional form of this expression resembles its counterpart for a usual two- 
dimensional medium. This might be explained from the dispersion relation Eq. ffTTl) 
where in the left-hand side we have something like the norm of a vector computed 
with respect to a metric with two spacelike components and a timelike one. Also, 
compared with a two dimensional photon gas, here the spectral density contains an 
extra factor, which is of the order of inverse characteristic length a -1 , leading to a 
huge enhancement of the density. The above relation is plotted in Fig. [3] which shows 
radical enhancement of the spectral density compared with a usual photon gas in 3 + 1 
Minkowskian space. 




Figure 3. Plot of the spectral density (in units of ^^4^ ) for X = 100 (upper) 
compared to the result of the usual spherical dispersion relation (lower). 



In spite of the above arguments, one may still be interested in studying the 
behaviour of the spectral density for situations where, say, X ~ 10 1 , corresponding 
to a microscale patterning length a ~ 10 _6 m. This is depicted in Fig. HJ Qualitatively, 
this looks like the previous case. 
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Figure 4. Spectral density (in units of 7r2 ^3a4 e 3 ) f° r X = 10. 



4. The partition function 

The thermodynamics of the above described photon gas can be derived from the 
partition function in the appropriate limit. Considering a grand canonical ensemble 
for the photon gas, we can write the grand partition function as follows 

pv r°° 

In Q = -— = - / p(u) ln(l - e-^)dco (16) 
Inserting Eq. (|12p into the above expression, we get 



+ 3 " 4 )-^/>}ln(l- e 



UJ 



)dw (17) 



where a ~ 0.67 is the root of 4 — 10x 2 + 3x 4 = 0. Now, using our previous argument 
on largeness of X = j3h~Q for a nanopatterned metamaterial, we can ignore the last two 
integrals in the above expression. This leaves us with 

in which 

f(x) = x (4 -10-^ + 3-^) ln(l-e-). 
In view of the largeness of X, this can be approximated as 

ln ^6^w/o e*-l (19) 

Inserting the values of the Bose-Einstein integrals above, this reduces to 

where C is the Riemann zeta function. 
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The internal energy is given by 



SV " '«3)-!k(5) + ^C(7)] (2!) 



3tt 2 c 3 /3 3 ^ V ^ 2 a 4 

which resembles the energy of a usual photon gas in two dimensions together with 
correction terms of the form X~ 2 and X~ 4 which are very small for a metamaterial 
medium with nanoscale patterning length. Comparing the last two equations, one 
reaches the following relation 

P = ~(\-0{X~ 2 )) (22) 

in which higher order correction terms are neglected. Other thermodynamical quantities 
such as entropy can also be obtained from the above partition function. The result of 
these calculations is consistent with the result obtained in the previous section according 
to which the photon gas under consideration behaves like an ordinary photon gas in 2 + 1 
dimensions with some small corrections. 

It would also be interesting to consider a metamaterial with e x = e y > and e z < 0. 
This can be modelled by interchanging e x and e z in Eqs. (|9]) and (flQj) . The relevant 
dispersion relation is then 

In this case, the wave-vector space configuration is again a hyperboloid, but with a 
different topology, namely a two sheeted hyperboloid depicted in Fig. [5j 




In this case, we are interested in the volume between the surfaces spanned by the 
relation Eq. (123 j) and the cutoff planes k z = ±k = j^\ / 2Q 2 — u 2 . Again, this volume 
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can decrease or increase with increasing values of u. The behaviour of the derivative of 
this volume with respect to frequency is shown in Fig. [6j 




Figure 6. The function v' = (in units of ^§j§-) in terms of lo (in units of f2) for 



the surface k 2 x + k% - ( % - 1) kl - ' J ~ - - : 



Repeating the same steps as in the previous sections, for u < Q we reach 

p(u) = v \-n 6 + n 4 u 2 + snV 

-3u 6 + {Q 5 - 3Q 3 u 2 )V2Q 2 - u 2 (24) 
This can be approximated as 

This shows a constant density of states corrected with a w 2 term. One may interpret 
this as the density of states of a photon gas in 1 + 1 dimensions for large values of Q. 

5. Discussion and conclusions 

We have considered a photonic gas with modified dispersion relation propagating in a 
nonmagnetic metamaterial medium with some negative components of the permittivity 
tensor. From the hyperbolic dispersion relations, we obtained the density of states, 
which in the case where the xx and yy components of the permittivity tensor are 
negative, shows some interesting features such as vanishing at a certain frequency. 
For the case where only the zz component is negative, the surface spanned by the 
dispersion relation in the wave-vector space is still hyperbolic but with a different 
topology and the density of states can decrease or increase with frequency. We studied 
the thermodynamics of the system and showed that it resembles a photonic gas in an 
ordinary one or two dimensional medium but with enhanced density of states. The 
results are of interest in studying the behaviour of physical systems in spacetimes 
with non- Lorentzian metric signatures. The relevant calculations are performed under 
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the assumption that the medium dielectric coefficients have a plasma-like frequency 
dependence but are independent of the wave- vector. It would be interesting to consider 
more general situations where these are functions of both the frequency and the wave- 
vector. 
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